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1. Introduction 

Precision computations beyond the leading order (LO) play a crucial role in the physical 
analyses performed at the LHC. Several new-physics signals, as well as many background 
processes, involve multi-particle final states. At LO, these processes exhibit a large scale 
uncertainty which can be softened by including the next-to-leading order (NLO) contri- 
butions, in which one-loop amplitudes have been a challenging component in processes of 
current importance [1,2]. 

Recent years have seen rapid progress in computing one- loop amplitudes, due largely 
to the use of on-shell methods grouped under the name of "generalized unitarity." One 
advantage of on-shell methods is that they enable the computation of loop amplitudes in 
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terms of tree- level amplitudes, which are relatively easier to generate either numerically or 
analytically in compact forms. 

Unitarity methods for one-loop amplitudes [3-18], recently reviewed in [19-21], depend 
on the knowledge of a small, canonical set of master integrals, in which any amplitude 
can be expanded linearly with coefficients that are rational functions of the kinematic 
invariants [22-37]. Generalized cuts of a loop integral are defined as the singularities in 
the limit of sets of virtual momenta going on shell: propagators are replaced by the delta 
functions of the on-shell conditions. For certain cuts, complex-valued momenta must be 
considered in order to solve all on-shell conditions. The coefficients of the master integrals 
are then obtained by matching generalized cuts of the loop amplitude and the master 
integrals. 

Unitarity methods work most beautifully when all internal particles are massless. If 
massive particles are involved, the same cuts can be used to produce coefficients of most of 
the master integrals without any new conceptual difficulties. More significantly, there are 
some additional master integrals whose cuts are more difficult to solve. In the case of QCD 
processes involving n massive quarks with masses mj = i v .. jTl , the additional master integrals 
are the tadpoles Ao(mi), and the "on-shell bubbles" Bo(m 2 ;mi,0) and Bq(0; m«, mi). The 
massless counterparts of these integrals vanish in dimensional regularization. These bubbles 
are called on-shell because they are in the momentum channel of a single on-shell external 
particle, either the new massive particle or any of the massless ones. There are several 
proposals for computing some of these coefficients [12, 21, 38-42] either numerically or 
analytically, with or without unitarity cuts. 

Here, in the context of renormalization, we will be concerned with the tadpole integral 
and the bubble Bo(m 2 ; m, 0), which has exactly the shape of the external leg correction 
diagram of the massive particle. Renormalized amplitudes also need the differentiated 
bubble, B' (m 2 ;m,0). The derivative is taken with respect to the first argument. 

External leg correction diagrams contain an internal on-shell propagator, which implies 
the singularity of their cuts, and hence the singularity of the cuts of the amplitude. This 
problem was raised and nicely addressed in [18,21], in the context of a numerical algorithm. 
The solution given in those papers was to omit the problematic contributions, taking care 
with the associated breaking of gauge invariance. 

Here, we present a different solution. In the spirit of the unitarity method, we want 
to keep the ingredients of the cut as complete amplitudes, without discarding any contri- 
butions. We must therefore also include the corresponding counterterms. To regularize 
the singularity arising from the cut of the external leg correction diagram, we perform an 
off-shell continuation of the cut momentum. The divergence is cancelled between the cut 
loop diagram and the cut counterterm diagram, so we retain only the finite terms in the 
expansion around the on-shell kinematics. The examples given here involve 4-dimensional 
cuts, although the formalism is equally valid in D dimensions. We are also able to sum 
over physical polarization states only, in the spinor-helicity formalism. 

The paper is organized as follows. In Section 2, we present the theoretical analysis. We 
define the off-shell continuation and show the cancellation of divergences for both bubble 
and tadpole coefficients from double and single cuts, respectively. In Section 3, we illustrate 
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the method with two examples computed with Feynman diagrams, namely H — >■ bb and 
qq tt. In Section 4, we illustrate the method in the spinor-helicity formalism for one of 
the partial amplitudes in it — > gg. For this calculation, we make use of compact analytic 
formulas for tree amplitudes. In Section 5, we summarize the algorithm. Appendix A lists 
the master integrals and some useful reduction formulas, and Appendix B gives identities 
and formulas for the massive spinor formalism used in Section 4. 

2. Description of the method 

The unitarity cut of the amplitude of a process with massive external particles should 
give information about the coefficients of the on-shell Green's function. Unfortunately the 
external leg correction diagrams are singular, because the propagator opposite the external 
leg carries the same momentum and is therefore also on shell. 

In order to apply finite unitarity cuts, we perform an off-shell continuation by doing 
a double momentum shift. We shift the massive external momentum k and one other 
external momentum r, to preserve momentum conservation. Assuming that both momenta 
are outgoing, the shift may be written as 

k^-k = k + £k, r->-r = r — (2.1) 

We choose the momentum k such that r ■ k = k 2 = and 2(k ■ k) ^ 0. The momentum k is 
on shell: k 2 = m 2 , while k is off shell. The momentum r is still on shell after the shift. With 
this shift, the propagator opposite the external leg diverges as l/£; thus the full amplitude 
diverges the same way. The cuts are calculated in terms of tree amplitudes. For the double 
cut, one has simply the three-point interaction Mr, now continued off shell, and an (n + 1)- 
point tree amplitude Ml, on shell but depending on the parameter £. In the case of the 
single cut we have a single (n + 2)-point on-shell tree amplitude Mt, which is ^-dependent. 
We expand Ml and Mt up to first order in £. The coefficients of this expansion could 
be obtained directly from a closed form for the full amplitude. Alternatively, separate 
recursive constructions could be used to generate them independently. 

The external leg correction diagrams are exactly cancelled by the corresponding coun- 
terterms, provided that the renormalization constants are defined in the on-shell scheme. 
The counterterm diagrams are constructed from the renormalization constants and the 
various tree-level off-shell currents. The renormalization constants are known in terms of 
the master integrals. The other ingredients needed for this step are the various tree- level 
off-shell currents. An expansion in £ is performed here as well, and the divergent part 
is simply discarded since it is guaranteed to cancel the cut loop diagram. The off-shell 
currents are gauge-dependent. In the sum of all parts, gauge invariance is restored by 
construction: coefficients of master integrals are gauge invariant, and we have only added 
zero in the form of the external leg correction plus its counterterm. 

After the cancellation of the divergent contributions, the on-shell limit is reached by 
setting £ = 0. Therefore, at every stage of the procedure, we systematically neglect terms 
of O(0- 
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Figure 1: Double cut of the external leg correction diagram, and the left and right tree-level 
amplitudes. The cut momenta are I and k — i. Color information is indicated by c and c' . The 
massive propagator giving the on-shell divergence is denoted by V . 

In the remainder of this section, we demonstrate the cancellation of divergences in the 
cut between loop Feynman diagrams and the counterterm. Our conventions for the master 
integrals and their cuts are given in Appendix |A|. 

2.1 Bubbles from double cut 

Consider the double cut of the external leg correction diagram for a massive fermion, as 
shown in Figure |l]. Let k denote the outgoing momentum of the external fermion, and £ 
and k — I the momenta of the cut gluon and cut fermion, respectively. 
The shifted massive propagator V is 

= ijjn + $ + = i(m + jt + $) 

(k + £k) 2 -m 2 £7 1 ' ' 

where 7 is defined by 

7 = 2k ■ k. (2.3) 
The tree-level amplitudes Ml and A4r depicted in Figure || read as follows: 

ML= (k + Ckf'-m 2 (%+^(™ + ^ + ^)A" Cc 
Mr = -ig T^ c , (u k u k+ f k _ e ) . (2.4) 

Here -4 c »c oxt denotes the remaining parts of the diagram, including the external legs on the 
left and all color-flow information. 1 Specifically, c" denotes the color of the propagator V . 

Note that the external massive spinor u k is not being shifted. Its shift is possible 
but unnecessary, since terms arising from the 0(£) contribution of u k drop out once the 
unrenormalized amplitude and the counterterm diagrams are summed together. We are 



1 Since it depends on k and f, -4 c " Coxt depends on £. 



- 4 - 



performing the complete momentum shift on Ml, supplemented by the specific instructions 
given above for how to continue the propagator V and the vertex Mr. 

Multiplying these expressions, we find that the double cut including the sum over 
polarization states, in Feynman gauge where Y1i £ ij. £ v = ~9^-> ls 



j dfi 2 k (2m 2 - £7) (u k Ac oxt ) + (u k (. (m + Jd + £$) Xc cxt ) + m£ (u k I A CCcxt 



2ig 2 C F 

The double-cut integration measure d/j, 2 t is defined in equation ( |A.2| ) , and the Dynkin 
index of the fundamental representation of SU(N) is 

2N 

It accounts for the internal color sum, leaving a Kronecker delta function allowing us to 
replace c" by c in A c " Ccxt . 

Using integral reduction (see ( |A.4j )), we get the bubble part of the diagram: 



2 C 

' Am 2 ( u k i CCext ) + 2£m ( u k I i CCcxt J B (m 2 + £7; m, 0) . 



167T 2 £ 7 L 

Finally, we expand around £ = 0, using 

B (m 2 + £ 7 ;m,0) = #0 (m 2 ; m, 0) + £jB' (m 2 ; m, 0) (2.5) 

and 

The result is that the bubble parts of the shifted diagram are given by 
Mb = -^^{4m 2 ^(u k A CCc jB' (m 2 ;m,0) 



167T 2 £7 

4m 2 (u k Ace**) + 4 ™ 2 £ («* ^c Ccxt ) + 2m £ ^ A 



CCcxt j 



Bo(m 2 ;m,0) L (2.7) 



When the amplitudes are written in the spinor-helicity formalism [43-47], the pro- 
cedure described above must be modified. Indeed, when using spinors, the completeness 
relation for polarization vectors is that of a light-like axial gauge rather than Feynman 
gauge: 

£\i £ v — —g^u H ^— • (2.8j 

A=± ' Q 

Here £ is the momentum of the gluon, and q is an arbitrary light-like "reference" momentum. 
The reference momentum (chosen independently for every gluon in the process) is needed 



to express polarization vectors in terms of spinors. Using eq. (2.8), we see that the double 
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Figure 2: Single cut of the external leg correction diagram. 



cut gets an extra contribution of the form 



ig 2 C F 



i 



CC cxt 



+ 



{u k f u k+fk _ e ) (u kHk _ £ I (m + Jt + £#) A Coxt ) 



The second of these terms vanishes in the sum over diagrams, 



(2. 



2 ^ H™ + $ + &) Ac cxt ) , (2.10) 

,4 

due to the Ward identity with the on-shell cut gluon. 
The remaining first term in (^.9|) is equal to 

«7 2 Cf /\ u fc I (m + # + ^ (m + # + £l) l CCoxt 
/ d ^2fe o ' ( 2>11 ) 



which gives a contribution to the bubble part of the amplitude, (2.5), that is 

g 2 C F UkH jijn + g Ac cxt 2 

16^ ~k B ° {m ;m ' 0) - (2J2) 

This new contribution does not affect the divergent part of the bubble coefficient. 

2.2 Tadpole from single cut 

Now we consider the single cut of the massive propagator of the external leg correction 
diagram. See Figure 2 The ingredients are nearly the same as in the double cut, 
except that the polarization sum for the gluon is replaced by the propagator —ig^y/l 2 . 
The tree-level amplitude Mt is 

2 C 

M T= ~ ^2^7 ^ k ^ U k+tk-l) {^k+fk-e 7m ( m + # + $) ^c'ccxt) • (2-13) 



2 We do not compute the single cut of the gluon propagator since the massless tadpole vanishes in 
dimensional regularization. 
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M 



ct 



Figure 3: Counterterm diagram. 



The single cut A. reads as follows: 



2g 2 C F 
£7 



(2m 2 - £7) {u k A CCc J + (u k t (m + $ + £$) A Coxt ) + m£ I A 



l cc cxt 



f 2 



The single-cut integration measure d/j, 1 t is defined in equation (|A.3[) . The integral reduc- 



tion relations ( A. 5 ) allow us to compute the tadpole portion of the external leg correction 
diagram. The result is 



9 2 C F 
16vr 2 



(2m 2 £7) (u k A CCc ^ + £m (u k % A 



"'•CCoxt 



^7(m 2 + £7) 



A (m) 



(2.14) 



Expanding in £ — >• 0, we get 
9 2 C F 



M A 



16tH 



A (m). (2.15) 



2.3 Cancellation against the counterterm 

Now we write the external leg counterterm, Ai ct , depicted in Figure ||, in terms of the 
master integrals, and we show that the divergent part is exactly cancelled by the divergence 
of the external leg correction diagram. 
The contribution of Ai is 



1 



M ct = - — (u k [j/iSZ^ + ilbZ^, - m5Z^ - m8Z m ) (# + $ + m) A c 



(2.16) 



In the on-shell scheme, the renormalization constants 5Z m and 5Z^ are expressed in terms 
of master integrals as 



sz m 

5Zjp 



9 2 C F 
16vr 2 
9 2 C F 
16vr 2 



A (m) 



m 
A (m) 



+ 2B (m 2 ;m,0) 
4m 2 B' Q (m 2 ;m,0) 



(2.17) 
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Therefore, the counterterm diagram is 



M ct = M rn + M i ', 



(2.18) 



where M m is given by 



M' 



9 2 C F 
16vr 2 



g 2 c F 

16vr 2 



+ 



+ 



+ 



2 4m 2 
— A (m) + -—B (m 2 ; m, 0) 

1 2m 

A (m) + 5 (m 2 ;m,0) 

7m 7 

2 4m 2 
— 4)M + -—B {m 2 ; m, 0) 

£7 £7 

2 4m 2 
-4)(m) + B (m 2 ;m,0) 

7 7 

1 . / v 2m „ , n 
A (m) + B (m 2 ;m,0) 

7m 7 



(2.19) 



while .M 1 ^ is given by 



9 2 C F 

167T 2 

9 2 C F 
16vr 2 



r A (m) + 4m 2 B' (m 2 ;m, 0) 



— 2 A (m) + 4m 2 B' (m 2 ;m, 0) 



m 



(2.20) 



by adding the following term to eq. ( 2.16| ): 



When the spinor-helicity formalism is used, the extra contribution ( |2.12| ) is accounted for 

+ ei + m)X Ccxt , 



M = - —u k 



(2.21) 



where 



SZ. 



, g 2 C F B (m 2 ;m,0) 



16vr 2 



k 



(2.22) 



Having expanded the counterterm around £ = at 0(£°), it is straightforward to check 
the cancellation of the terms proportional to and the ones which depend on k once 
■M-A, Mb, and M ct are added together. The sum is identically zero, since in the on-shell 
scheme the external leg correction diagram is exactly compensated by the external leg 
counterterm. The actual contribution to the tadpole and bubble coefficients comes from 
the other diagrams in the full amplitude. They are finite in £, thus no k dependence arises 
in the £ — > limit. 
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b(k 2 ,c 2 ) 




K fc 3,C3) 



b(k 2 ,c 2 ) 



H(fci) 




6(^3, C3) 



H(fei) 



b(k 2 ,C 2 ) 



b(k 3 ,c 3 ) 



Figure 4: Virtual corrections to the process fl3.1|). 



3. Examples 

3.1 The process H ^ bb 

As a first application, we focus on a simple process, the Higgs decay into a bottom-anti- 
bottom pair: 

H(h) -> 6(fe 2) pa) 6(fc 3 , c 3 ). (3.1) 

We apply the procedure described in Section |2| to compute the coefficients of the on-shell 
bubble Bo(m^;mb,0), of B' Q (ml;mb,0) and of the tadpole Ao(mb) entering the NLO QCD 
amplitude M-H^bb [48-51]. We regularize the tree-level amplitudes using the following 
shift: 

k 2 ^k 2 = k 2 + k 3 ^k 3 = k 3 - £k. (3.2) 

The massless momentum k is such that 

72 = 2(k 2 -k)^0, 73 = 2(k 3 -k)^0, (3.3) 

while k 2 and fc 3 are on-shell: k 2 = k 3 = m%. Here we use a slightly different shift since 
both the shifted momenta are off shell. In this example, the modification of the method is 
irrelevant and streamlines the computation. 



The coefficient b 
the double cut, 



B {ml+~i2£,;m b fi) 



of Bo(ml + 72^; m;,,0) can be computed by means of 



d H,k 2 M H^bgb\ K ^ 



i,k 2 -£,£,k 3 ) M gb ^ b {e,k 2 -£,k 2 ), 

where the tree-level amplitudes on the r.h.s. are related to the processes 

H(ki) -> b{k 2 ,c 2 ) g{k A ,A 4 )b(k 3 ,c 3 ), g(h,Ai)b(k 2 ,c 2 ) -> b(k 3 ,c 3 ) 
respectively. Up to 0(£), the relevant part of the double cut is given by 



m (u k2 (.Vk 3 ) + (h ■ i) (u k2 v kz 
(k 3 + £) 2 - ml 



+ 



8g 2 em b 5 c c l 

3M\yS w 

{uk 2 1 %. + $) v k3 ) +m (u k2 (/ + % + £%)v k3 ) - - mf ) (u k2 v ks ) 



272^ 



+ 



(3.4) 



.(3.5) 
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Using the results ( |A.4j ) , we get the coefficient 



g 2 em b 5%i 



Am 2 b - M H ^ 



m b (2m b + 372C) r \ <t(- % m I /q o\ 

(m 2 b + l2t)l2t i mt ^^ ) + ^"^^)J C (3 ' 6) 

The parameter e is the charge of the positron, while Mw and s w are the mass of the W- 
boson and the sine of the weak mixing angle respectively. 

The coefficient b Bo ^ m 2 _^ 3 . mb<0 - ) of Bo(ml — £73;?rif>,0) can be obtained analogously 
starting from the double cut 

\k 3 M H^bl = J d ^2h M H ^ bgb -(h,k 2 ,£,k-£) M gb ^ b (£,k-t,h), (3.7) 
and using the relations (|A.4j) . M. gb ^ b is the tree- level amplitude of the process 

g(k 1 ,A 1 )b(k 2 ,c 2 ) -> b(k 3 ,c 3 ). 

The result is 



_ g 2 em b 5%i 
b B ( m l- 13 c,m b ,o) ~ 24tt 2 M w s w 



2(2m 2 b - M 2 H ) _ 2m 2 + l3 ^ 

-Ml 73? [Uk2Vk3) 



. m b {2,^-2m 2 ) f , % . , 1 

+ 73 g( m 2 _ 73 g) ^ I"* 2 * V k- A ) + wt 6 (n fc2 ^jj >• (3.8) 

The coefficients of Bo(ml;m b ,0) and of B' (ml;m b ,0) can be obtained by expanding 

6 BoK+ 72?;mb , o)Bo(ml + 72?; m&, 0) + b Bo{m 2_^. mbfi) B (ml - 73?; m 6 , 0) (3.9) 
around ? = and neglecting terms of 0{^). 

The tadpole coefficient a^ r mb \ can be computed from the sum of two single cut dia- 
grams, 

A ljt2 M H ^ bb + A 1 _~ k ,M H ^ bb = J d/ii, M Hb ^ bbb (k u £,k 2 ,£,k 3 ), (3.10) 

provided that the loop momentum i is chosen appropriately. 3 The tree-level amplitude on 
the r.h.s. is related to the process 

H(k x )b(k 2 ,c) -> b(k 3 ,c 3 )b(k 4 ,c)b(k 5 ,c 5 ). (3.11) 



In particular in each diagram of Figure ^ the loop momentum has to be fixed such that the only internal 
propagator appearing are £ 2 , {£ + k^) 2 — m 2 , and (£ — k 3 ) 2 — m' 2 .. 
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6 

b(fc3,C3) 



6(^3, C3) 



Figure 5: External legs counterterm diagrams of the process ( p.l[ ). 

Divergent diagrams with zero-momentum internal gluons do not contribute, owing to the 
color flow of the process fl3.11| ). At the part of the single cut ( |3.10| ) which is relevant 

for the computation of the tadpole coefficient reads as follows: 

Ag 2 em b 5 c c 2 i f J (u k2 1 {kf 2 + £%) v ki ) + m b (u k2 / v ks ) 
' dfj, 1Q ' 



3M w s w J [ (£ + k 2 + £A;) 2 72 £ 



(l+k 3 - ik) 2 73^ 



+ ■■■}. (3.12) 



The coefficient a j 4 ( m(j ) is obtained using the relations ( |A.5| ) and expanding the following 
expression in £: 



g 2 em b 5^i 



2Air 2 M w s v 
1 



(u k2 %v k3 ) + m 2 (u k2 v ks )] 



(m( + 72O72C 



+ (^-^)(n fc2 , fc3 ) >. (3.13) 



The results ( |OD , (|T|), and (p^) , have b een checked performing Passarino-Veltman 



decomposition of the one loop amplitude obtained from the Feynman diagrams in Figure ^j. 
This procedure has been performed with the help of FeynArts [52], FormCalc [53,54] and 
FeynCalc [55]. 

We have verified the cancellation of the 0{£,~ l ) terms against the contribution of the 
two counterterm diagrams depicted in Figure [5], along the lines of Section 2^ . In the £—>•() 



limit, the terms containing k drop out. The bubble and the tadpole coefficients are thus 
independent of this unphysical momentum. 

3.2 The process qq — > ti 

A less trivial example involves the process of top-anti-top production via quark-anti-quark 
annihilation: 



q(k 1 ,cl)q(k 2 ,c 2 ) -> t(k 3 , c 3 ) i(fc 4 , c 4 ), (qj^t). 



(3.14) 




q 9 





Figure 6: NLO QCD corrections to the process ( 3.14 ). The ghosts are labelled by rj while q' 
u, d, c, s, t, b. 



This is one among the channels entering the hadronic production of a top-anti-top pair 
whose NLO QCD contributions have been computed [56-59] and implemented in MCFM [60] 
and in MCONLO [61]. We have computed the coefficient of BQ(m^;mt,0), -B^mJ 2 ; mt, 0) and 
of the tadpole ^4o( m t)- The tree- level amplitudes have been regularized using the shift 



^3 ^3 — ^3 + ki — ?> fej = ki — £k. 



(3.15) 



We have used the double cuts 



\~k 3 M qq->tt = J d ^2,h M qq->tgt(h,k 2 , h - £,£, fc 4 ) M gt ->t(£, h ~ £, h), (3.16) 

\k 4 M qq^tt = / dfi 2 k 4 M qq ^ tg i(k 1 ,k 2 , fa,£, h ~ t) Mgt->i(ti ^4 - h), (3.17) 



the single cut 



*V. M "+* + \-i* M "+e = f^ M q ^ ta (k M k,iM), (3.18) 

and the relations ( |A.4j ) and ( |A.5| ). The coefficients can be obtained performing an expansion 
around £ = at along the lines of Section 3.1. The actual expression of the coeffi- 

cients is lengthy and is omitted. We have checked our result against the Passarino-Veltman 
based computation of the Feynman diagrams depicted in Figure ||. In the renormalized 
amplitude, we have checked the cancellation of the 0(£ -1 ) contributions and of the terms 
containing the massless momentum k. 
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4. On-shell bubbles in the spinor-helicity formalism 



In this section we illustrate our method in the spinor-helicity formalism. We compute 
the bubble coefficient of an amplitude from the usual cut integral using double cut spinor 
integration [62-67]. 

As discussed in Section |2[ the spinor-helicity formalism requires a gauge choice of a 
null "reference" momentum for each gluon. On-shell quantities are independent of this 
choice, but it plays a role in our method because of the off-shell continuation. Of course, 
the total bubble coefficient will be gauge- invariant. However, the independence of gauge 
choice will not be apparent until all counterterm diagrams are considered, including the 
internal ones. 

The cut integral involves the off-shell continuation of a three-point amplitude which 
is a varying function of q, the reference momentum for the cut gluon. Furthermore, the 



off-shell current A CCaxt (defined in Section 2A) depends on the reference momenta chosen 
for the external gluons. 

In practice, one makes a convenient choice of reference momenta. It is important to 
be consistent in this choice throughout the entire calculation. The gauge choice of the 
reference momentum q of the cut gluon cancels out immediately between the cut and its 
counterterm. In the following example, we keep q general to show this cancellation, while 
making convenient choices for the reference momenta of external gluons. 

The spinor-helicity formalism can accommodate massive particles as follows [45,68]. 
For a particle of mass m and momentum p, the massive spinors are given by 

[p v t]\ \p v rj) 

[7] p p J ir] p"J 

where rj is a fixed null momentum, and we define the null vector 

p =p - - rj. (4.3) 

2p ■ Tj 

For an antiparticle, the mass m should be replaced by — m everywhere. 

It follows that massive spinor products obey the following (anti) symmetry relations: 

(ij) = - 0'»> . M = - . <v] = b*> ■ ( 4 - 4 ) 

See Appendix ^ for further details. 

As an example, we now focus on the process 

t(ki,a,hi) t(k2,c 2 ,h 3 ) -)• g(k 3 ,c 3 ,h 3 ) g(k4,C4,h4), (4.5) 

where hi (q) is the helicity (color) of the i th particle. We will compute the coefficient of 
Bo(m^;rnt,0) of the leading-color part of the all- minus helicity amplitude. Contributions 
to the on-shell bubble coefficients come from the unrenormalized amplitude M-j^ in the 
k\ and k\ channels and from the counterterm diagrams. 
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Figure 7: fc^-channel cut of the right and left primitive amplitudes. 
4.1 Contributions of the unrenormalized amplitude 

At one loop and in D = 4 dimensions, the color decomposition of the unrenormalized 
amplitude reads as follows [69,70]: 



tt-^gg U n \ 



N{T«T*)%AW(l*\2!»,i**,jt> 

(M)e{(3,4),(4,3)} 



+ tr (T C3 T C4 ) S 1 9 2 ^ , 3^ 3 , 4^ 4 



(4.6) 



where N is the number of colors. The one-loop partial amplitudes are the leading-color 
ones, while A^l is subleading. 4 

The external leg corrections of the anti-top are leading color, thus we focus on the 
leading-color partial amplitudes. We compute the bubble coefficient of A^\ (if , 2~[ , 3~ , 4~) . 

The coefficient of A^\ (if, 27", 4~, 3~) can be obtained from the previous one via the sub- 
stitution 3f>4. The leading-color partial amplitude can be written in terms of primitive 
amplitudes as [69] 

41 (If >2 t -,3-,4-) = A L (If ,2T,3- 4-) - ±A R (if, 2 t ~, 3,", 4") 

+ ^E4' [1/2l (l^2 t -,3 9 ",4-), (4.7) 

where the sum runs over the quark species. The primitive amplitudes contributing to the 
coefficient we are interested in are the left amplitude A L and the right amplitude ^4^. 
They can be obtained starting from the color-ordered Feynman rules [71] and selecting the 
diagrams without a closed fermion loop, where the fermion flow passes to the left (A L ) or 
the right (A R ) of the loop. 



4 The factor of 2 is related to the different normalization of the color matrices used through the paper 
and the ones in the color-ordered Feynman rules as in [71]. 
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4.1.1 The cut in the k\ channel 

We perform the shift setting k = k 4 , i.e. 



ki -»• ki = h + £k A , fc 4 -> k 4 = (1 +£)k 4 , 7 = (4|1|4]. 



(4.8) 



Note that these shifts have the same sign, since k\ is incoming while k 4 is outgoing, in 
contrast to our all-outgoing convention of Section |||. 

The double cut of the right amplitude A R , Figure 0, reads as follows: 



T R 



2,-fci 



(4.9) 



The integrand I R is written in terms of color-ordered tree-level helicity amplitudes as 



i R = Y, A(i h -{t + h)T h \-rMAUt + h)f,2i,?,-A- 9 J7 



(4.10) 



ht,h g =± 



The necessary helicity amplitudes are collected in Appendix [B,2 



After the spinor integration and the £ expansion, the bubble coefficient ^ ( m 2 + ^ 7 . mt ) 
can be written in a relatively compact form: 



1 2ml (34) [12] mj (34) (m t [2q] (qi) [41] + [14) <g[l|2] [4g]) 



B {m 2 t +il\mu0) £ r^i (4|1|4] 2 

m t 



+ 



[34] (4|l|4] 2 (g|l|g] 



[34] (4|1|4]< 



mt (41] [42] (34) + m|(32](41] + m t (4|1|2] [13) 



2 [12] (34) (4|1|4] 



(4.11) 



The dependence on q comes from the off-shell three-point function. See equations (|B 
and ( B.lOj ). This dependence will drop out in the sum with the counterterm diagram. 
The double cut of the left amplitude A L reads as follows (cf . Figure 0) : 



r 

2,-fci ' 



2,-fci 

with the integrand I L defined as 



i L = y A(l t - > -4 A ^-(^+fc^)^' it )A((^+fc 1 )J*,4^2 t -,3J,4J 



(4.12) 



(4.13) 



h t ,h g =± 



By considering the properties of the color-ordered Feynman rules, the qgqgg amplitude 
can be expressed as a linear combination of qqggg amplitudes, as depicted in Figure || 
Therefore I L can be written as I L = I R + 1 , where I is given by 



I l = A^-ie + h)^,-^) L((£ + fe 1 )^,2r,3-,^,4-) 



+ A((l + k 1 )^,2t,e h g \3-Ag 



(4.14) 
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I 6 



- k-, / 4 





Figure 8: Pictorial representation of the relation between the qgqgg and the qqggg color-ordered 
tree-level amplitudes. 



It turns out that the bubble coefficient from I 1 vanishes. Therefore, the bubble coefficient 
of the left amplitude, ^ ( m 2 + ^. mi \ is equal to the one of the right amplitude: 



h L - h R 

°B (m? +£ 7 ;m t ,0) U B (m t 2 +£ 7 ;m t ,0) m 



(4.15) 



The final contribution to the bubble coefficient of the leading color part of the ampli- 
tude (4.6) is obtained using the results ( 4.11| ) and ( 4.1 5| ) and expanding Bo(m 2 + £7; rnt-, 0) 
around £ = 0. The outcome is 



2m\ (34) [12] 2 



b B (m?Hr,rn u o) B (m 2 ;m t ,0) 



[34] <4|1[4] 



B (m t ;m t ,0) 



+ (3^4) . 



(4.16) 



4.1.2 The cut in the k 2 channel 

In this case, the symmetry of the amplitude allows us to write down the result as a simple 
relabeling of the contribution from the k\ channel. It suffices to exchange 1 O 2, 3 <4 4, 
and mt —> —mt- The result is 



2.9 



i g 2 C F 

1 8vr 2 



(T C3 T C % 



b R 



B (m 2 ;m t ,0) 



2m? (34) [12] 2 



B (m t ;m t ,0) 



+ (3 O 4) 



[34] <4|2|4] 

4.2 Contributions of the counterterm diagrams 

The counterterm diagrams can be color decomposed as follows: 

■2 



(4.17) 



M 



ct 
tt-tgg 



(r*r*)ZA« (if ,2f ,3f,4f) + (T^)^ (lf,2f,4j,3j) 



. (4.18) 



The color-ordered counterterm diagram A ct can be computed using the color-ordered Feyn- 
man rules after performing the shift fl4.8| ). We will focus on A ct (if , 2f , 3~, 4~) , since 
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A ct (l f ,2 t ,4 g ,3 g ) can be obtained via the substitution 3f>4. The corresponding color 
ordered counterterm diagrams are depicted in Figure ||. 

The diagram D a in Figure |] is related to the tree-level amplitudes where lj is con- 
tinued off-shell. These amplitudes can be derived by a recursion relation as described in 



Appendix B.2. Using a decomposition similar to the one in eq. ( |2.18| ), the counterterm 
diagram D a can be written as D a = A m + A^ + A k , whose bubble terms can be read 
from eq. (|2l9|) , ( g3g ), and (|2T2l|) . Substituting the expansion of J from (|B.17| ) into these 
expressions, we find 



^ <4|1|4] 2 [34] 



8vr 2 



B (m 2 ;m t ,0) 



m 2 t [l\3\2] 
(4|1|4] [34] 5 



+ B (m 2 ;m t ,0) 



mt (34) ([14] (4|1|4] [32] + m t [14) [43] [42]) 



(4|1|4] 2 [34] 2 



B (m 2 ;m t ,0) 



A k 



9 2 C F 
8vr 2 



m t [1|4|3|2] 
2 <4[1|4] [34] 5 

2m t 3 (34) [12] ) 
<4|1|4] [34] I 



9 2 C F 
8tt 2 



B (m 2 ;m t ,0)- 



m\ (34) (m t [q2] [14] (4q) + <<j|l|2] [14) [q4]) m t [12] (34) 



(4|l|4] 2 ( (? |l|g] [34] 



(4|1|4] [34] 



We have kept the full (/-dependence above, in order to demonstrate that it cancels exactly 
between the bubble computed from the cut and the counterterm. In practice, we might set 
q = throughout, in which case the g-dependent terms simply vanish. 

The external leg counterterm diagram in Figure ^ can be obtained from D a by 
performing the substitutions 1 2, 3f>4, and mt — > —mt. 

Finally, there are contributions to the coefficient of Bq (m 2 ; mt , 0) coming from all 
other self-energy counterterms. In this case, there is only one, namely diagram D c . The 
counterterm separates the scattering process into two off-shell currents. For the example 
at hand, the currents are simply the cubic interactions, and the Feynman rules give the 
bubble coefficient 



9 2 C F 



m t 



8vr 2 (4|1|4] 2 [34] 2 



2m? - (4|1|4] [14) [34] (32] + [13] (43) [42] 



+ 2mt(4|l|4]([14)[42]-[13](32] 



(4.19) 



The bubble contribution of the counterterm diagrams, eq. (P^) , is then 



2g 2 



(T C3 T^)H (D a + D b + D c ) + (3^4) 



(4.20) 
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Figure 9: Counterterm diagrams contributing to the coefficient of the on-shell bubble 
B (m^;m t ;0). 

4.3 Total bubble coefficient 

Here, we assemble all the contributions to the coefficient of -E>o(?ri£ ; mt, 0). The bubble 
coefficient coming from the cuts and counterterms is obtained summing the contribu- 
tions ( [4.16| ), ( |4.17| ), and ( |4.20| ), After the cancellation of divergent and gauge-dependent 
pieces between the cut-diagram and counterterm contributions, and then recovering the 
on-shell limit, we find the following result: 



To compute the bubble coefficient for the entire amplitude, one would need the con- 
tributions from all remaining helicity amplitudes. It is easy to translate the results above 
to any other helicity amplitude whose gluons have equal helicities, either by overall parity 
conjugation of the amplitude, or by conjugating the individual massive spinors. The con- 
ventions for massive spinors (see Appendix [B]) allow their individual conjugation except in 
spinor strings where they appear on both ends; these strings can be split using Schouten 
identities. Thus the remaining nontrivial computation for the leading-color part would be 
any single configuration with opposite helicities for the two gluons. 

There are also coefficients multiplying the differentiated bubble, B' (m^;mt,0), arising 
from the various counterterms via 5Z^. They are proportional to the tree- level amplitude. 
The counterterm for external wavefunction renormalization cancels exactly with the cut 
diagram in this coefficient, as shown in Section [|, but other counterterms need to be treated 
separately. 

5. Summary 

We close this article with a brief review of our procedure. 

1. For each massive external fermion, with momentum denoted by k, choose the shift 
direction k. Define the off-shell continuation k = k + ^k, where £ parametrizes 
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the shift, and choose another external leg to absorb the shift in order to conserve 
momentum, as in equation fl2.1|) . 



2. Compute the tree- level inputs for the cuts in the shifted momentum channels. The 
cuts may diverge as l/£. The expressions are needed through 0(£ ' ). The expansion 
might be done by a contour integral, or by generating the 0(£ -1 ) and terms 
separately. 

3. The cut integrals have now been regulated, and the coefficients of master integrals Aq 
and Bq can be found in the usual way. Keep only the 0(£°) term. The 1/^-divergent 
term is ignored, unless it is desired for a consistency check. It will be proportional to 
the tree amplitude, as shown in equations ( |2.7|) and ( |2.15 ). Because of the off-shell 



continuation, the cut depends on the gauge choice made for the cut gluon. 

4. Compute the counterterms from tree-level currents. The counterterms for external 
leg corrections are computed at The divergent part cancels the divergence of 

the loop exactly and may likewise be ignored. Similarly, the coefficient of B' will 
cancel and may be discarded. In the spinor-helicity method, one must include the 
term A4k from equation ( |2.21| ), making the same gauge choice for the cut gluon as 
in the cut integral. 

All other counterterms for internal legs and vertices are finite and can be computed 
from unshifted expressions. Considered separately, the counterterms depend on the 
gauge choices made for external legs. Gauge invariance is restored in the sum. It is the 
internal counterterms that provide the surviving contributions to the differentiated 
bubble B' , via the renormalization constant SZ^,. 
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A. Master integrals, double and single cuts 

The master integrals that pertain to the wavefunction renormalization diagram are the 
scalar bubble, 
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and the tadpole with massive propagator, 



C27 r//) ; n 

iit 2 J (k — i) 2 — m 2 



Mn, )= ^— J d D i 



These integrals are dimensionally regularized by taking the dimensionality D = 4 — 2e. The 
double cut of the bubble and the single cut of the tadpole are obtained by substituting the 
propagators with the corresponding delta functions: 



A 2 , k B (k 2 ; m, 0) ee (2?r ^ 2 J d D £ 5(i 2 ) 5 ((k - if - m 2 ) , 

A 1:k A (m) = {2lT ^l ° j d D i 5 ((k - if - m 2 ) . (A.l) 

VK J 

The rule for computing the double cut from the left and right tree amplitudes is 

A 2 , k M = J dfi 2 , k Ml(£)Mr(£), dii 2 , k = -T^u d D i 5(£ 2 ) 5 ((k - if - m 2 ) .(A.2) 

The rule for computing the single cut from tree amplitudes is 

A 1>k M = f cKfc Mt(£), d^ k = dPi 5 ((k - if - m 2 ) . (A.3) 

For convenience, we list the reduction of several double-cut integrals, 

J dn2,k= ~ A 2 , k B (k 2 ;m,0), 

J d^ k (R ■ i) = - -±J^—^1{R . k) A 2 , k B (k 2 ; m, 0), 
f , (R-1) i k 2 (p-R)-(k-p)(k-R) A „ /i2 . , A 

and single-cut integrals, 

7 ^ k Ji + pf = ~ 32^ A i.^o(m). (A.5) 

Here, i? and p are any 4-vectors. 

B. Massive spinors 

We use the formalism of massive spinors developed by Kleiss and Stirling [45]. 

For assistance in calculations, we have used the specific representation of massive 
spinors in terms of an arbitrary null "reference" momentum rj, which is the same for all 
particles, as follows [68]. For an on-shell particle of mass m and momentum p, define 



b m 



2 



p =p - 77, (B.l) 

2p ■ 7] 
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which is another null vector. 

Then the massive spinors are defined as 



{# + m)\rj\ (fl + m) \r]} 

\p) = r b ) , \P] = i $ v , (B.2) 
[p v rj\ \p v T]) 

whose conjugates are 

[rj\(p + m) (r/| (jf + m) 

(P = — f yi—i iP\ = —/ jA~ • ( B -3) 

For an antiparticle, the mass m should be replaced by — m everywhere. These formulas 
are smooth in the massless limit. 

We have also used the Mathematica package S@M [72] to help with spinor manipulations 
and numerical evaluation. 

B.l Some identities for massive spinors 

Among the identities obeyed by the massive spinors, we find the following particularly 
useful in our computations. 

The contractions of spinors of the same chirality are antisymmetric as usual, while the 
mixed spinor product is symmetric: 

(ij) = -(ji), [ij] = -[ji], (ij] = \ji)- (B.4) 

The spinor representation of a momentum vector can be replaced by the massive spinors 
using 

f = \p) [p\ + \p] (p\ ~ 771 f° r a particle, (B.5) 
f = \p) \p\ + \p] (p\ + rn for an antiparticle. (B.6) 

The slash is implicit inside spinor products. Note that = 2p. L -pj. 

Finally, the Schouten identity is still valid for any four spinors of the same type, 

(ij) (kl) + (ik) (Ij) + (il) (jk) = 0, [ij] [kl] + [ik] [Ij] + [il] [jk] = 0, (B.7) 

while of course the three-spinor uncontracted versions remain valid for any massless spinors 
a, 6, c, 

|a) (be) + \b) (ca) + |c) (ab) = 0, [a] [be] + |6] [ca] + \c) [ab] = 0. (B.8) 
B.2 Tree- level formulas 

Here we list the tree amplitudes needed for the cut in the example of Section ||. They 
are adapted from [73]. 5 The three- and five-point amplitudes, which are ingredients in the 
cut computation, are given in a form that allows straightforward parity conjugation of the 
fermions. 

The convention we follow here is that gluon momenta are directed outward, while 
quark and anti-quark momenta are directed inward. The mass of the quark is to. 



^Equivalent formulas were derived earlier in [74-76], also by on-shell recursion relations. 
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Three-point amplitudes. The relevant three-point amplitudes, valid also in the off- 
shell continuation, are 

A(lr,2-,3J) = [l|( Mia±mi )|2], (B.9) 

A(l.,2^3j)=[l|(Ml±jaM)|2]. (B.10) 

The quantity inside parentheses is the polarization vector, which depends on the null 
reference momentum q. On-shell, the amplitude is independent of q. 

Four-point amplitudes. The tree-level amplitude corresponding to the one-loop am- 
plitude of interest is 



m (34) [12] 
(2p 2 -p 3 )[MY 



A(ir,2-,3-,4") = -fi^ . (B.ll) 



Five-point amplitudes. The amplitudes entering the double cuts ( [4 . 9[ ) and ( [4.12 ) are 

m (3|p 45 ll|5)((14] [32] -(13] [4 2] ) 
(2p 1 -p 5 )(2p 2 -p 3 )[34} 2 [45] 

[5|l|p 34 |2[3) ([25] <1|P34|2|3) - [2|p 3 4|d 5 i |3) (15]) 



2" 3;,4J,5J) = 1 n ',2r„, ~ . (B.12) 



A(\g > 2„ , 3„ , 4„ , 5 



^(1^ j 2 , 3 q , 4+, 5 ff 



(2pi.p 5 )(2j>2-jJ3)[34] [45] (5| P 34|2|3) 

m (43) 3 ((l|p 34 |5) (3M ~ (lNj> <5|p34|2]) 
(pi+p 2 ) 4 (45)(35)(5|p 34 |2|3) 

m (35) 3 ((l|p34|5) (3|p 4 5|2] - (l|p 45 |3) (5|p 34 |2]) 



(B.13) 



(pi+p 2 ) 4 (34) (45) (5| P 34|2|3) 
[4|1|5> [4|2|3) 



m (14] (53) [42] 



(2pi-p 5 )(2p 2 -p 3 )[34] [45] (5|p 34 |2|3) 
(15> [42] [4|2|3) - <14] <32] [4]1[5> V (B.14) 



AiXij ■* 2^ , 3 fl , 4 g , hg 



3|2|p 45 |l|5) ((13] (5|l|p 45 |2] + (l|p 45 |d 32 |5) [32]) 

(B.15) 



(2pi-p 5 )(2p2-P3)[34] [45] (5|l|p45|3> 
m(45) 3 ((l|p 34 |5) (3|p4s|2] - (l|p4s|3) (5|p 34 |2]) 



(Pi+p 2 ) 4 (34)(35)(5|l|p 45 |3) 

where pij = pi +Pj and dij = Pi — pj. Similar amplitudes with the opposite helicity choice 
for l<j are given by exchanging angle brackets |1) , (lj with square brackets |1] , [1|. (The 
quark 2 q can be conjugated similarly, but we do not need that in the present paper.) This 
type of conjugation is valid as long as the derivation of the tree amplitudes did not apply 
any identities such as the Schouten identities, (B.7), which does not respect conjugation of 
a single massive spinor if it is contracted with another. 

Off-shell four-point currents. For the counterterm, we use a more general formula for 
the tree- level four-point current. This expression is valid when the particle lq is continued 
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off shell and its external spinor is stripped off. Because it is off shell, it depends on the 
reference momenta q 3 and 54 of gluons 3 and 4, respectively. The full expression is 



1 



[9 3 3] [? 4 4] 



1 



<3|2|3] 



[lg 4 ] <4|1|? 3 ] (32] - [14} [ 94 |2|3> [g 3 2] 



m [lg 4 ] (43) [q 3 2] - m [14) [q A q 3 ] (32] 



+ — ( [lg 4 ] (42] [q 3 4] + [14) [g 4 2] [g 3 4] 



+ [lg 3 ] (32] [ 94 3] + [13) [q 3 2] [q A 3} + 



[l|d34|2] fe^] 



(B.16) 



For the example discussed in Section ^j, we need the current for the process where only 
k\ is off shell, and we choose #3 = & 4 , q 4 = k 3 throughout. This current, with the spinor [1| 
stripped off, is given by 



J 



p u |2|3) [32] + m |3] (43) [42] + (m 2 - p\) |3] [23) Pu 1 2] 



2p 2 - P3 [34]' 



2 [34] 



2 • 



(B.17) 
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